Weakly bound states of neutrons in gravitational fields 
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Abstract 

In this paper a quantum-mechanical behaviour of neutrons in gravitational fields is considered. 
A first estimation is made using the semiclassical approximation, neglecting General Relativity, 
magnetic and rotation effects, for neutrons in weakly bound states in the weak gravitational field 
of the Earth. This result was generalized for a case, in which the Randall - Sundrum correction to 
Newton's gravitational law on the small scales was applied. Application of the results to Neutron 
Star physics is considered and further possible perspectives are discussed. 
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I. INTRODUCTION 



The idea concerning the potential of neutron storage by using of UCN (ultra cold neu- 
trons), probably was first suggested in 1959 by Ya.B. Zeldovich [1]. In his work a simple 
estimate concerning the possibility of UCN conservation in the container was devised. It 
was shown that a UCN (with wave length larger than 500 A and effective temperature less 
than 10~ 3 K) should be totally reflected from the reservoir walls (which were made from 
either carbon or beryllium material) when the neutrons are at velocities around 5m /s and 
less. 

More rigorous results, obtained later, completely confirmed these estimations. In general, 
these results can be obtained in the framework of pure quantum - mechanical or optical 
approaches. In the quantum mechanical approach the potential of the reflecting surface 
is constructed on the basis of the average procedure of the pointlike Fermi quasipotential, 
given in the form: 

u(r) = • bS(f— fj) 

m 

where b - is a neutron wave length of coherent scattering on the nuclei of the considered 
surface. In this approach, by strictly following the framework of quantum mechanics it 
was shown that the depth of neutron penetration in the surface material are significantly 
less, than the neutron wave length (2- 4 ] and the neutron scattering on the surface can be 
considered as elastic. An alternative approach 5], was based on the optical analogy of light 
scattering on metals, such that all inelastic processes in the reflection can be effectively 
described by the imaginary part of the refraction coefficient. In this case the work of 1] 
serves a necessary theoretical premise for the starting point of all further experimental works 
based on the UCN technique. One of the first experimental results was obtained in 1968 by 
the JINR (Joint Institute of Nuclear Research, Dubna) {f| . The importance of this type of 
work is extremely high, because it opens a new and unique possibility to carry out highly 
precise experiments, such as accurate measurement of neutron lifetime or its electrical dipole 
moment [7] and etc. 

The first proposal using experimental measurements with UCN in the gravitation field of 
the Earth was done by [8] . Therefore, there appears to be sound knowledge base from which 
to proceed given the achievements to date in both the theoretical and experimental fields. 
Certainly, the many recent works devoted to this research subject indicate that the solution 



of the UCN storage problem can give us an unique capability to carry out a very precise 
experiments in the field of neutron physics and neutron interferometry in gravitational fields 



In this paper, we will begin by reviewing some pertinent recent experimental results 
and measurements of UCN weakly bound quantum states in the weak gravitational field of 
the Earth 
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15| . After that, we will present a simple sketch of some of the other possible 



applications of these results to astrophysics and gravity. A more detailed description of some 



of other applications can be found in 



16]. 



II. ENERGY OF WEAKLY BOUND STATES OF NEUTRONS IN THE EARTH'S 
GRAVITATIONAL FIELD. 



Theoretical consideration of a neutron's energy spectrum and its wave functions in the 
Earth's gravitational field, as cited in the 15), can be found in a number of works, including, 



for example 



ljl9 



19l | . In this section, by using pure methodic reasoning, we present a simple 



and transparent derivation of the main relationships of these results. A description of weakly 
bound states in the gravitational field of the Earth follows in the Schrodinger Equation 
(SE) solution using a simple spherically symmetric case with gravitational potential energy, 
defined as: 



6V{r) 



-jMm(l-±\ r> R 



(1) 



oo, 



< r < Rn 



where M and m are the Earth and neutron masses at rest respectively and r - is a distance 
from the Earth in cm up to the point where the neutron mass m is placed. In this case SE 
can be written in the following way: 




E - 5V(r) 



tp = 



(2) 



where ip is a neutron wave function and E is its corresponding energy. Substituting, in the 
normal way, ip = into the equation (2) we have: 



d 2 f 2m 
— - H 

dr 2 h 2 



E - 5V(r) 



f = o. 



(3) 
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Using a simple approximation for the gravitational potential energy in the vicinity of the 
Earth surface, we can write: 



SV(r) ->■ SV{x) = 7 



Mm 

Rq 







X I 


X \ 











(4) 



where r = R + x, ^Cl and Rq is an average radius of the Earth. In this approximation 
SE can be written as: 



dx 2 



+ (a - (5x)f = 0, 



(5) 



where we introduced a and (3 as parameters, defined by: 

2m ^ ( m\ M 



(6) 



By substituting a new variable z = a — fix, we can finally write SE in the convenient form 
for the analytic solution: 



d 2 f 



r/ = 



(7) 



dz 2 (i/3) 2 ' 

For this ordinary differential equation of the second order, as suggested by the theory of 
Bessel functions [20] , there exists a regular solution which is presented in the following form: 

2zi' 



f(z) = C N y/Z 



MW} +J -H 3/3 



(8) 



where is a normalization coefficient of the wave function. From this last expression we 
can extract the energy spectrum of the neutrons weakly bound states by using a boundary 
condition at the point x = 0, which simply leads to the equation: 



, , 2at\ , /2a? 







(9) 



If you assume that ( n — > ( 



3 

2a 2 

3/3 



is such a number, then it satisfies the equation: 



J,[( n )+J_,(( r 



(10) 



which leads us to: 



E„ 



9m 



Co(Cn 



(11) 



TABLE I. Comparison of our results with theoretical and experimental results, obtained in [131 ]. 
Energy values are given in peV units and the Cn value is dimensionless and E n -th. is calculated 
using expression (17) 



n Cn E n -th E n ,\13\ -exp. E n , our 

1 2.3834 1.3767 1.4 1.4054 

2 5.5105 2.42191 2.5 2.4573 

3 8.6474 3.27356 3.3 3.3180 

4 11.7868 4.0255 4.1 4.0794 

5 14.9272 4.71261 - 4.7751 




FIG. 1. Cn roots distribution for the Ji (Cn) + ^_i(Cn) = equation. 



where we introduced cq = ^-^(Hg) 2 j • By now using approximate values for the above 

defined constants, such as: g « 9.80655 ms" 2 ; R « 6.371 • 10 6 m; mc 2 « 939.565330 MeV 
and He ~ 197.327053 MeV- fm, we can calculate the first bound states for the E n values and 
compare them with other results, including experimental measurements. These results are 
presented in Table 1, where the Co value is set to Co ~ 0.757325 • 10~ 12 eV. 10mm In Fig.l 
an approximate distribution of the Cn roots is presented. 
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III. ASYMPTOTIC PROPERTIES 



Given the solution obtained above it appears that it would be interesting to analyze its 
asymptotic properties. To do this we are able to use the properities of Bessel functions {20) ]: 



Uv) 



2 v 

— COS 

Try J 



y- tt(^ 



(12) 



by introducing y — > we can obtain: 



f{y) ->■ ^y/n^y e sin 



7T 



(13) 



from which we arrive at: 



/(*) 



z 4 sin 



2ql 
~3£~ 



(14) 



For the zero roots of this function we can immediately obtain from the last expression that: 

(15) 



2a§ 



7T 

3/3 + 4 



727T 



and finally, using our notations for variable z we are led to: 



E n = mgx + 



9/ " fc \2V 

ling) I I // — 



(16) 



Naturally, one might find the last expression to be quite interesting. This is true because 
the 1-st term of the expression can be precisely viewed as the classical part of an energy 
particle in a homogeneous gravitational field while the second term is a quantum mechanical 
contribution. In the case where h — > we have arrived at a pure classical result and in the 
case where H ^ 0, but x = we have exactly obtained results matching the energy spectrum 
from |l3| : 



En 



1 



2 

1\ 3 

/# " _ 4 



(17) 



IV. THE SEARCH FOR HIGH DIMENSION CONTRIBUTIONS. 



As it was pointed out in |13[ , the possibility of making observations concerning the UCN 
bound states in the Earth's weak gravitational field can provide new insight regarding the 
verification of known interactions involving small distances. 



An interesting view on the gravitation interaction, theoretically can lead us to the idea, 



elegantly declared in 
lar, it follows from 



21] , regarding its universality and mult i- dimension nature. In particu- 



22] that a high dimensional correction using the framework of the Randall 
- Sundrum (RS) theory can be applied to Newton's gravitational law at small distances. To 
obtain a better understnding of this premise it is useful to evaluate such theoretical predic- 
tions from the experimental point of view. An interesting proposal related to the search of 
;he RS correction to Newton's gravitational law on small scales was recently suggested by 



23]. 

In this section our goal is to theoretically estimate the possible contribution the RS 
correction might have on Newton's gravitational law at small distances on the UCN weakly 
bound states in the Earth's gravitational field and to calculate an upper limit of the Irs 
(Randall - Sundrum) parameter: 



V(r) 



- r fMm 



Mi + *]-i(i + ^f 



oo, 



r>Rn 



< r < R 



o 



Using an approach simular to the one applied in the previous section we find that in a 
spherically symmetric case that the equation (3) can be rewritten by using the ^ < 1 
approximation as: 

d 2 f 2m 



+ 



E - nuj \ 1 + 3 • % |.r 



f = 



(19) 



dx 2 h 2 

Introducing in a similar way as before we inject a new additional constant: Ars = 7^-, 
therefore the last equation (19) can be rewritten as: 

d 2 f 



dx 2 



+ (a- (3x)f = 



(20) 



with redefinition of the a and (3 values from equation (5) it follows: 

2 



a = a 



2mE 



/3 = 2| - I <?(1 + 3A 



2 

RS 



(21) 



As it was shown before the solution of the equation (20) can be written in the form: 



f(z) = c N Vz 



Ji 



2zi 
3/3 



+ J_ 



255 
3/3 



(22) 



where z = a — (3x. Using a simple connection between (3, and $ values we arrive at: 

P = 0(l + 3A 2 J] (23) 



further we can decribe the contribution of the A RS term regarding of weakly bound states 
of the UCN in the Earth's gravitational field using the boundary condition at x = 0: 

.'.r:V':f? ; i m 



3 



3£ / M 3/3 



From the previous equation by substituting a new variable ( = ^ ) the final results can 
be presented as: 

E n = h^(9Kn) 2 (l + 3AL) 2 ) 3 = c (l + 3A^ 3 (I (25) 

where ( n are a root of the equation (24). From this expression it follows that A RS is equal 
to: 

from which we can estimate, that 

l RS < 3.9 • 10 3 m (27) 



V. QUANTUM MECHANICS APPLIED TO THE NS SURFACE. 

Currently, one actual problem of great import in the area of astrophysics is the detection 
the Black Holes (BH) which could serve as an experimental proof of Einstein's gravitational 
theory. Another very important problem is the development of different precise methods for 
the measurement of the physical properties of distant astrophysical objects, such as: quasars, 
pulsars, gamma ray bursts, neutron stars (NS) and etc. In all these cases it would be very 
useful to provide such experimental investigations by different and completely independent 
physical methods. These experiments if successful would immediately raise the validity of 
the experimental data and provide clarification concerning the physical picture of the objects 
under investigation. Specifically, the greatest impact would occur in BH and NS physics. 
Within NS physics there may be a very interesting connection regarding the interplay of 
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general relativity (GR) and the quantum mechanical effects in the NS interior. The existence 



of the strange stars (SS) 24| underline that NS is a research area that is a very interesting and 



important for our further understanding of the behaviour of matter at extreme conditions. 



In this section we begin with a simple theoretical estimation because of the [llNl5| approach, 
for the quantum effects on the NS surface. Introducing new notations: M^s - NS mass, 
Me - Earth mass, we can derive using the framework of the nonrelativistic approach of 
gravitational acceleration on the NS surface which is equal to: 



Mns\ ( Re 



2 



where Re and Rns are the Earth's and the NS radii respectively and rj is defined as 77 



;/ 1 1 , N , . If the density of neutron matter on the NS surface is p = mn p , then we 

can simply estimate an average size of the cell, which would contain one neutron as: 

1 

/ m\ 3 , h , . 

a = - > A c = 29 

VP/ m^oc 

where m r o - is a 7r° meson mass, p - is the neutron matter density on the NS surface, n p - its 
concentration, A c is the Compton wavelength and a - the size of the cell. In cases in which 
the cell size is larger than the 7r° meson Compton wavelength, in that it is approximately 
equivalent to the short range of nuclear forces, we can conclude that neutrons on the NS 
surface can be considered in the first step approximation as a free gas in the external grav- 
itational field of the NS. This is because the average distance between them is larger than 
the radius of the nuclear interaction. 

In our consideration, we suggest, that due to the Pauli exclusion principle there should 
exist, close to the NS surface, an analogy of the Fermi like surface, which can play the 
role of a mirror as is the case of the experiments conducted with the UCN in the Earth's 



gravitational field 



ll|4l5| 



In this case, for rough estimation, we can simply follow the expressions derived before for 
the Earth's gravitational field with small modification to obtain the final result for the NS. 
The final expression, after such a modification of expressions (11), would appear as follows: 



1 
3 

2 „/c 2 



Ei N 5) = [^(anghfj ■ r/i E^ifi (30) 

Finally, a simple example to illustrate the calculations described earlier related to the NS 
with given parameters is provided, the chosen parameters follow: M NS = 2 • M Sun ~ 4 ■ 10 30 



TABLE II. Comparison of the energy of neutron's bound states (in eV units) on the NS and Earth's 
surface 



n E n on NS E n , on Earth (th.) 

1 5.68 -KT 5 1.40 • 1(T 12 

2 9.92 • 1(T 5 2.46 • KT 12 

3 1.34 • 1(T 4 3.32 • 1(T 12 

4 1.65 -10" 4 4.08 • 1(T 12 

5 1.93 -10" 4 4.77 • 10~ 12 



kg, R NS fa 10 4 m and M E fa 5.98 -10 24 kg, R E fa 6.37- 10 6 m. The results of such calculations 
are presented in Table 2. For these numerical data rj fa 2.7 ■ 10~ n . From the data in Table 
2 we can conclude, that perhaps it is possible to discuss the possibility of NS spectroscopy, 
where the transition wave lengths will be around \t r ~ ~ (0.124 — 1.24) cm, and 

the according transition energy AE tr ~ (10~ 5 — 10~ 4 ) eV. Details of this spectroscopic 
information are directly connected to the NS physical parameters, such as the M^g and 
Rns values. The probability of their spectroscopic transitions n — >■ k can be described as 
an overlap of the corresponding wave functions as described by the simple relation: 

u n ^k = (^J ip n (r)ip k (r)drj (31) 

VI. CONCLUSIONS 



The present research, which used the works of [Ill4l5j] as a foundation devised a simple 
manner in which to consider the energy of UCN weakly bound quantum states in the Earth's 
gravitational field. Comparison of our results with the results of the above cited papers is 
presented in the Table 1. 

The obtained results were generalized using the case of Newton's gravitation law correc- 



tion by applying the RS theory 22j in the framework of a nonrelativistic approach. It was 
determined that in the further development of measurement techniques of the UCN in the 
Earth's gravitational field, we can more precisely estimate the upper limit of the contribution 
of the RS correction to Newton's gravitational law at a small distances. The contribution 
to the UCN energy spectrum from the RS correction to the formation of their bound states 
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was derived and the upper limit was estimated as Irs < 3.9 • 10 3 m. The ramifications of this 
are that the experimental methodology probably is not sensitive enough to check for a RS 
type correction to Newton's gravitational law at small distances. 

The second step involved the estimation of neutrons bound state formation in the external 
gravitational field of a NS in comparison to the Earth's gravitational field neglecting GR, 
magnetic and rotation effects. In this case the results obtained are presented in Table 2, 
which presents some sample NS parameters and compares them with theoretical results 
obtained in the Earth gravitational field. 

It is obvious, that when considering the neutron bound state formation on the NS surface 
it is necessary to apply a complete GR evaluation of the problem. The first and transparent 
approach in this direction is a reconsideration of the SE solution for the formation of a 
neutron bound state in the form: 

g ij {v i V j - r*.W)v + - E/ e //(r))v = (32) 

where g v is a metric of the gravitational source (for example NS) and U e ff(r) is its effective 
external gravitational field. Here i,j, k are spatial indexes. More mathematically rigorous 
considerations would need to use the Dirac wave equation in the gravitational field of the NS. 
Of course this approach would need to take into account magnetic fields and rotation effects, 
but is beyond the scope of this work in which the goal was to present some preliminary 
theoretical estimations. Application of such an approach can be very useful because of 
the potential of considering NS radial oscillation 25] which can lead to changes in the 
gravitational field on the NS surface. Also, there may be some observable effects related 
to neutron spectroscopy. In conclusion, we want to mention the paper 26], where the 
UCN storage problem was considered using the previously mentioned magnetic mirror in 
the presence of a gravitation field. 
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